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Abstract. The structure theory for the quadratic algebra generated by first and second 
order constants of the motion for 2D second order superintegrable systems with nondegene- 
rate (3-parameter) and or 2-parameter potentials is well understood, but the results for the 
strictly 1-parameter case have been incomplete. Here we work out this structure theory 
and prove that the quadratic algebra generated by first and second order constants of the 
motion for systems with 4 second order constants of the motion must close at order three 
with the functional relationship between the 4 generators of order four. We also show that 
every 1-parameter superintegrable system is Stackel equivalent to a system on a constant 
curvature space. 
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1 Introduction 

A classical second order superintegrable system H = £V ■ Q^PiPj+Vfa) ° n an n-dimensional local 
Riemannian manifold is one that admits 2n — 1 functionally independent symmetries £fc(x • p), 
k = 1, . . . , 2n — 1 with C\ = H, that are at most second order polynomials in the momenta p\. 
(Further, at least one Ch = aa^(x)piPj + TUh(x) with h > 1 must be exactly second order.) 
That is, {7i,Ck} = where 



{f,g} = ^(dxj fd Pj g - d Pj fd Xj g) 



3=1 

is the Poisson bracket for functions /(x, p), g(x, p), Here 2n — 1 is the maximum possible 
number of such symmetries, For the case n = 2 The structure of the Poisson algebra generated 
by the symmetries has been the subject of great current interest. For potentials depending 
non-trivially on 2 or 3 parameters, see [I], for a precise definition, it has been shown that the 
algebra is finite-dimensional and closes at order six in the momenta. All such algebras have 
been classified, as have been all spaces and potentials that give rise to them [2JG]. Similarly all 
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degenerate 1-parameter potential systems are known via case-by case classification as well as 
the associated Poisson algebras, [3]. The number of true 1-parameter, not just a restriction of 
a 3-parameter potential, systems is 15 (6 in complex flat space, 3 on the complex 2-sphere and 
one for each of the 4 Darboux spaces). Under the Stackel transform that maps superintegrable 
systems into equivalent systems on other manifolds, these divide into 6 equivalence classes. 
However information about the structures of the corresponding algebras are known only by 
a case by case listing and the mechanisms by which they close have never been worked out. 
Some results for 1-parameter potentials were reported in pQ but although the results are correct 
the 1-parameter proofs are incomplete. Here we work out the structure theory and prove that 
the quadratic algebra generated by first and second order constants of the motion for systems 
with 4 second order constants of the motion must close at order three and must contain a Killing 
vector. Furthermore we show that there must be a polynomial relation among the symmetries 
at order four. It is important to develop methods for understanding these structures that can 
be extended to structures for n > 2 where the analysis becomes more complicated, and this 
approach should point the way. 

We treat only classical superintegrable systems here, though the corresponding (virtually 
identical) results for the quantum systems follow easily [5]. In both the classical and quantum 
cases the symmetry algebras and their representations have independent interest [6j [HJ O 
[lot fill fl2l [TBI Q3J [15]. Also, there are deep connections with special functions and orthogonal 
polynomials, in particular Wilson polynomials \16\ I17|. 

As an example, consider the classical Hamiltonian on the two sphere 

n = j? + Ji + Ji + ^, 

s 3 

where the Ji are defined by <J3=sip S2 —s 2 p Sl and cyclic permutation of indices, and s\+S2+s\=l. 
If we seek all first and second order constants of the motion for this classical Hamiltonian, we 
find three possibilities in addition to H itself viz. 

_03_ , 2 2^ 4 a- a 3SlS2 

2*2 



= J? + ^2-{ 1 + 4-sl), A 2 = JiJ 2 -^^, X = J 3 . 



'3 °3 

The set of 4 symmetries X 2 , 7i, A\ and A 2 is linearly independent, but functionally dependent 
via the fourth order identity 

Ai{n - a 1 - x 2 ) - a\ - 1 {X 2 + n) + 1 = 0, 

They satisfies the Poisson algebra relations 

{X,A 1 } = -2A 2 , {X,A 2 } = -H + X 2 + 2A 1 , {A 1 ,A 2 } = -X{2A 1 + a 3 ), 

so the algebra closes at order three. We will show that this structure is typical for all 1-parameter 
potentials that are not just restrictions of 3-parameter potentials, that is there is always a Killing 
vector (a first order constant of the motion), the algebra always closes at order three, and there 
is always a fourth order relation between the 4 generators. 

The situation changes drastically for the two sphere Hamiltonian with nondegenerate poten- 
tial 

2 2 2' 

s l s 2 s 3 

where s 2 + s\ + s 2 = 1. The classical system has a basis of symmetries 

^2 ^2 ^2 ^2 ^2 ^2 

C\=Ji + CL 2 — + a3~2 ' £>2 = J 2 + a 3~ + °l — 2? £3 = J% + a l~2 + a 2~2 ' 

s 2 s 3 s 3 s 1 s 1 s 2 
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where Ti = £1 + £ 2 + £3 + 0,1 + 02 + 03 and the Ji are defined by J% = sip S2 — S2p Sl and cyclic 
permutation of indices. The classical structure relations are 

{£1 , 11} = 8d (H + a x + a 2 + a 3 ) - 8£ 2 - 16£i£ 2 

- 16a 2 £ 2 + 16a 3 (H + ai + a 2 + a 3 - £1 - £ 2 ), 
{£ 2 , 7£} = -8C 2 (H + a x + a 2 + a 3 ) + 8£| + 16£i£ 2 

+ 16ai£i - 16a 3 (H + ai + a 2 + a 3 - £1 - £ 2 ), 

with {£x, £ 2 } = 1Z and 

ft 2 - 16£i£ 2 (H + ai + a 2 + a 3 ) + 16£ 2 £ 2 + 16£i£^ + 16ai£? + 16a 2 £| 

+ 16a 3 (H + ai + a 2 + a 3 ) 2 - 32a 3 (W + ai + a 2 + a 3 )(£i + £2) + 16a 3 £ 2 

+ 32a 3 £i£ 2 + 16a 3 £ 2 — 64aia 2 a 3 = 0. 

Now there is no longer a first order symmetry but 3 second order symmetries. The algebra 
generated by these symmetries and their commutators now closes at order 6, [18]. The com- 
mutator 1Z cannot be expressible as a polynomial in the generators, but 1Z 2 and commutators 
of 1Z with a generator can be so expressed. This 3-parameter system is called nondegenerate. 
Note that our 1-parameter potential is a restriction of the nondegenerate potential, but that the 
structure of the symmetry algebra has changed drastically. 

On the other hand, the system with the 1-parameter potential 

a a a 
v = — + — + —' 

s l s 2 s 3 

i.e., the restriction of the nondegenerate potential to the case a = a\ = a 2 = a 3 has a symmetry 
algebra that is exactly the restriction of the algebra for the nondegenerate case. Further, any 
2-parameter potentials obtained by restricting the nondegenerate potential can be shown not to 
introduce symmetries in addition to those obtained by obvious restriction from the symmetry 
algebra for the nondegenerate case. We will show that these examples are typical for 2D second 
order superintegrable systems and will clarify the possible structures for the symmetry algebras 
in the degenerate cases. 



2 Background 

Before proceeding to the study of superintegrable systems with potential, we review some basic 
facts about second order symmetries of the underlying 2D complex Riemannian spaces. It is 
always possible to find a local coordinate system (x,y) = (xi,x 2 ) defined in a neighborhood of 
(0, 0) on the manifold such that the metric takes the form 

ds 2 = X(x, y){dx 2 + dy 2 ) = A dz dz, z = x + iy, z = x — iy, 

and the Hamiltonian is Ti = (p 2 + p 2 ) /A + V(x, y) = TCq + V, where V is the potential function. 
We can consider a second order symmetry (constant of the motion) as a quadratic form £ = 
2 

a l: >(x, y)piPj + W(x, y), = a? 1 , that is in involution with the Hamiltonian 7i: {7i, £} = 0. 

'..7 = 1 

2 

A second order Killing tensor £ = ^(x, y)piPj is a symmetry of the free Hamiltonian Ho: 
{TCo,Cq} = 0. The Killing tensor conditions are 

d? = -%a"-%a» i = l,2; 

2a^+a^ = -^ 1 -^ 2 , i,j = 1,2, i^j. (1) 
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From these conditions we easily obtain the requirements 



2a 



12 



,11 



2a 



12 



,11 



,22 



From the integrability conditions for these last equations we see that 



Aa 



12 



0. 



A(a 



ii 



,22 



)=o, 



a = dl + dl 



In order for a form C = Co + W to be a symmetry of the system TC = TCq + V it is necessary 
and sufficient that Co be a Killing tensor and that W satisfy the equation 

{H ,W} + {V,C } = 0. 

The conditions for this are 

2 



Wi 



1,2, 



where Wi = d Xi W, Vj = d Xj V. Necessary and sufficient that these last two equations can be 
solved is the Bertrand-Darboux condition 



22 



V n )a 12 + V 12 (a n 



,22 



) 



A 



V X + 



A 



For a second order superintegrable system we demand that there is a Hamiltonian and 2 
other second order symmetries: TC, C±, C 2 such that the Killing tensor parts of the 3 symmetries 
are functionally independent quadratic forms. By a change of basis if necessary, we can always 
assume that C\ is in Liouville form, so that the coordinates x,y associated with TC, C\ are 
separable. Thus, we can choose our orthogonal coordinates x, y = x%, x 2 such that the quadratic 

1. In this system we have A12 = 0. A second symmetry 



22 



,11 



11 



,22 



1/A, a = 0, which clearly always satisfies 



form in C\ satisfies a 12 = 0, a 
is defined by the Hamiltonian itself: a 

equations ([I]). Due to functional independence, for the third symmetry C 2 we must have a 12 7^ 
and it is on this third symmetry that we will focus our attention in the following. Now the 
integrability conditions can be rewritten as 



A 



12 



0. 



A = A 



22 



A 



11 



3Ai4i + 3X 2 A 2 - (An + A\ - A 22 -A 2 2 )\ = 0, 



(3) 



where A = In a , the subscripts denote differentiation and A satisfies An + A 22 + Af + A\ = 0. 
Equivalently, 

A12 



0, a\j + all = 0, 



a 12 (An - A 22 ) + 3Ai< - 3A 2 a^ 2 + « - a%)\ = 



.12 



12 



.12 



.12^ 



(4) 



In this second form a fundamental duality becomes evident |19t Q] (with a typo in the second 
reference): If X(x,y), a 12 (x,y) satisfy (jl]) then 



X(x,y) = a 1 



12 ( x + iy —ix — y 

7T' V2 



a (x,y) = X 



x + iy —ix — y 
v 7 ^ ' V2 



also satisfy these conditions. Thus, the roles of metric and symmetry can be interchanged, and 
a second interchange returns the system to its original state. 

Another key equation is the integrability condition derived from consideration of A12 = 0: 



5L( 1 )A 1 - 5L< 2 )A 2 + " 4 2) + 3^i£ (1) - 3A 2 L^)X = 0, 



(2) 



(5) 
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where L^> = An 2 — A\ 2 A\, = A\ 22 — A\ 2 A 2 . We will derive this in detail in Section[5j The 
dual version of condition ([SD is the integrability condition. 



bK^af - bK^af + (K^ - K™ + 3 Pl K^ - 3p 2 K^)a 12 = 



(2) 



(6) 



where p = lnX and = p 2 22 - 2pnp 2 - P22P2 + pjp2, = -pui + 2/9 2 2Pi + P11P1 - p\p\- 
Note that = K^> = is the condition that A is a constant curvature space metric. Indeed, 
this is exactly the necessary and sufficient condition that A(lnA)/A = c where c is a constant. 

Koenigs [19] employed condition © to show that the only spaces admitting at least 6 li- 
nearly independent constants of the motion were constant curvature spaces. Indeed in that case 
there are 3 functionally independent symmetries a 12 (x,y), hence 3 equations ([6]). That is only 
possible if the coefficients of a\ 2 , a2> 2 and a 12 vanish identically. Hence = = and 
A is a constant curvature metric. Koenigs did not make use of condition ([3]), but from our point 
of view this condition is more fundamental. 

Using our special coordinates and the Killing equations ([1]) we can write the Bertrand- 
Darboux equations (|2|) in the form: 



12 



A2 
A 



Ai 
A 



V 2 , 



V: 



22 



Vi 



11 



2Ai 
"X" 



+ 3Aj 



V\ 



2A2 
A 



+ 3^2 



Vo 



(7) 



In [T] we have shown that this system admits a nondegenerate i.e. 3-parameter (the maximum 
possible) potential V(x, y) if and only if the potential is the general solution of the canonical 
system 



V 12 = A"(x,y)Vi + B L \x,y)V 2 , V 22 - V n = A 22 (x, y)^ + B 22 (x, y)V 2 



(8) 



whose integrability conditions are satisfied identically. Thus this system will admit a 4-dimen- 
sional solution space, with one dimension corresponding to the trivial addition of an arbitrary 
constant. Each solution is uniquely determined at a point (xo,yo) by prescribing the values V, 
Vi-> V 2 , Vn. In our special coordinates we have 



A 12 = 



A2 
A ' 



B 



12 



A_i 
A 



^, A 



22 



B 



22 



(9) 



We will say little about 2-parameter potentials other than pointing out, as we already showed 
in p] that they are just restrictions of 3-parameter potentials. Their canonical equations take 
the form 



V12 
Vu 



A"(s,y)Fi+B"(x,y)^ 2 , 
A u (x,y)V 1 + B u (x,y)V 2 . 



V 22 = A 22 (x,y)V 1 +B 22 (x,y)V 2 



By relabeling coordinates, if necessary, we can always assume that the canonical equations 
for 1-parameter potentials take the form 

V 1 = B\x,y)V 2 , V 22 -V 11 = B 22 (x,y)V 2 , V 12 = A n B 12 (x, y)V 2 , 



where the integrability conditions for these equations are satisfied identically. This system will 
admit a 2-dimensional solution space, with one dimension corresponding to the trivial addition 
of an arbitrary constant. Each solution is uniquely determined at a point (xq, yo) by prescribing 
the values V, V 2 . In our special coordinates we have 



B 



12 



A2 R i Ai 



B 



22 



2^ + 3A l 
A 



B 



(10) 
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3 The Stackel transform 

The importance of the Stackel transform in superintegrability theory is based on the following 
observation. Suppose we have a superintegrable system 

\{x,y) 

in local orthogonal coordinates, with fc-parameter potential V(x, y), < k < 3 and suppose 
U(x,y) is a particular choice of this potential for fixed parameters, nonzero in an open set. 
Then the transformed system 

n = i±Pi + v(x,y), ~X = \U, V = y 

is also superintegrable. Indeed, let S = Yl a^PiPj + W = So + Whea second order symmetry 
of TC and Sjj = ]P a^piPj + Wjj = Sq + Wjj be the special case of this that is in involution with 
(Pi + U. Then it is straightforward to verify that 

Wjj 1 

S = S --^H + V H (11) 

is the corresponding symmetry of TC. Since one can always add a constant to a potential, it 
follows that 1/U defines an inverse Stackel transform of TC to TC. See [20} |2T| [22] for many 
examples of this transform. We say that two superintegrable systems are Stackel equivalent if 
one can be obtained from the other by a Stackel transform. Note from (jlip that the off-diagonal 
elements a 12 = a 21 of a symmetry remain invariant under the Stackel transform. 

If V(x, y) is a nondegenerate potential, i.e. the general (4-dimensional) solution of canoni- 
cal equations (jSJ) and U(x,y) is a particular solution of these equations, then V(x, y) is also 
a nondegenerate potential satisfying the canonical equations 

V22 = V n + A 22 Vt + B 22 V 2 , Via = A 12 ^ + B 12 V 2 , 

where 

i 12 = ^12_^2 ^22 = A 22 + 2 ^1 B 12 =B 12 -^, B 22 = B 22 - 2^. 

U : U' U' u 

Similarly, if V(x, y) is a 1-parameter potential satisfying canonical equations 

Vi = B 1 (x,y)V 2 , V 22 -V 11 = B 22 (x,y)V 2 , V 12 = B 12 (x, y)V 2 , 

and U(x,y) is a particular nonzero instance of this potential then V(x,y) is also a 1-parameter 
potential satisfying the canonical equations 

V 1 = B 1 (x,y)V 2 , V 22 -V n = B 22 (x,y)V 2 , V l2 = B 12 (x, y)V 2 , 

where 

B 12 = B 12 -2B^, B 22 = B 22 + 2((B 1 ) 2 - l)^, B 1 = B 1 . 

Note that the function B 1 remains invariant under a Stackel transform. 

Now we return to study of a general superintegrable system with /c-parameter potential, and 
viewed in our special coordinate system. Then A is the metric and V is the general solution 
of the Bertrand-Darboux equations ([7]). If the integrability conditions for these equations are 
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satisfied identically, then the solution space is 4-dimensional. Otherwise the dimensionality 
is less. The metric A must satisfy the fundamental integrability conditions © that depend 
only on a 12 , invariant under the Stackel transform. Now if U is a particular solution of the 
Bertrand-Darboux equations ([7]) then it defines a Stackel transform to a new Riemannian space 
with metric fi = XU. Since a 12 is invariant under the transform the fundamental integrability 
conditions for fj, are that same as for A: 

/ii2 = 0, H22- fin = 3/ii^i-3^ 2 ^2 + {A u +Al-A&-A%)ii. (12) 

Note that these two equations appear identical. However they have different interpretations. 
The fixed metric A satisfies ([3]) and is a special solution of (fl2j) . Here \i designates a k + 1- 
dimensional family of solutions, of which A is a particular special case. It follows that A satisfies 
the integrability conditions for this system. There is an isomorphism between the solutions U 
of d?]) and the solutions f/, = XU of ([3]). The \i are precisely the metrics of all systems that can 
be obtained from the space with metric A via a Stackel transform. 

4 Review of known results 

For the most symmetric case, the potential is nondegenerate so k = 3. Then the solution space 
of (|12p is 4-dimensional. In [2] we used this fact to derive all metrics and symmetries that 
correspond to superintegrable systems with nondegenerate potential. The possible symmetries 
A = In a 12 that can appear are precisely the solutions of the Liouville equation Ayi = Ce A 
where C is a constant. If C = then the system is equivalent via a Stackel transform to 
a superintegrable system in flat space. If C ^ then the system is Stackel equivalent to the 
complex 2-sphere. As is easy to verify, these systems are precisely the solutions of the system 
of equations 

L {1) = A U2 - A X2 A X = 0, L (2) = A 122 - ^12^2 = 0. 

An amazing fact is that these systems are exactly the same as those derived by Koenigs in his 
classification of all 2D spaces admitting at least 3 functionally independent second order Killing 
tensors. 

A referee has called our attention to two recent and very interesting papers by Tsiga- 
nov |23|. [24] . He assumes that a superintegrable system admits an orthogonal separation of 
variables in some coordinate system, so that there is a related Stackel matrix. Under this as- 
sumption one can construct the action angle variables as explicit integrals. Using these Tsiganov 
shows a rough duality between superintegrability and functional addition theorems, which al- 
lows one to find a generating function for constants of the motion, including those higher than 
second order. The Euler addition theorem for elliptic functions leads to the construction of all of 
the 2D superintegrable systems with nondegenerate potential! In this sense, these potentials are 
implicit in the work of Euler. This is clearly a powerful method for constructing superintegrable 
systems. It doesn't appear to yield any proof of completeness of the results or any classification 
of all possible spaces that admit superintegrability and in distinction to say [2] it requires the 
assumption of separation of variables. 

In [l] we showed that all 2-parameter potentials were restrictions of nondegenerate potentials. 
Further, assuming the correctness of Koenigs' results we carried out a case by case analysis over 
many years to find all superintegrable systems with 1-parameter potentials. We found that 
all of these were restrictions of nondegenerate potentials again. However, in some cases the 
restricted potential admits a Killing vector, so that the structure of the associated quadratic 
algebra changes. These results were not based on a theoretical structure analysis. Our attempt 
at a structure analysis, contained in [1] was incomplete and there were gaps in the proof, though 
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the results are correct. Thus there is reason to use our Stackel transform approach to revisit 
this issue for 1-parameter potentials. 

Finally the O-parameter case deserves some attention, although it was already treated by 
Koenigs. If we consider a O-parameter potential system as one in which it is not possible to 
admit a nonconstant potential, then it follows from our results and those of Koenigs that no 
such system exists. Koenigs' proof used complex variable techniques, very different from the 
methods used here. 



5 3-parameter and 2-parameter potentials 

Suppose we have a 2D second order superintegrable system with zero (or constant) potential, 
i.e., 3 functionally independent second order Killing tensors. Under what conditions does there 
exist a superintegrable system with nondegenerate potential V such that the potential-free parts 
of the symmetries agree with the given Killing tensors? To answer this we choose the special 
coordinates such that the integrability conditions for the zero potential case are: 

Ai 2 = 0, A = A 22 - An - 3AiAi + 3A 2 ^ 2 - (A u + A\ - A 22 -Af)X = 0, 
A = A n + A 22 + A\ + A\ = 0, 

where A = In a 12 . Necessary and sufficient conditions that this system admits a nondegenerate 
potential V with canonical equations 

V22 = Vu + A 22 V X + B 22 V 2 , V12 = A 12 V X + B 12 V 2 , (13) 

where A 13 , B %3 are given by Q, are that the integrability conditions for equations f)13j) are 
satisfied identically. These conditions are: 

T (1) = 2B\ 2 - Bf - 2A 12 - Af = 0, (14) 
T (2) = 2B\ 2 A 22 - A 22 Bf - A 22 A 12 - 2A 12 B\ 2 - A 22 + A% + 2A l2 A 12 

+ B 12 A 22 - B 22 A 12 - B 22 Af - A\ 2 = 0, 
y(3) = -B l2 Af + 2A l 2 2 B 12 + B 22 B\ 2 - B 22 Bf - 2B 12 B\ 2 - A 22 B\ 2 

+ A X2 Bf + B\l - B\\ - B\\ = 0. 

Substituting expressions ([9]) into (j!4[) we find that = identically. To understand the 
remaining conditions we use A = 0, A = to simplify the equations. We solve for Am and A222 
from Ai = 0, A 2 = 0, respectively, and substitute these expressions in and T^ 2 \ respectively. 
Then we solve for An — A22 from A = and substitute this result into and T^ 2 \ Then we 
find 

T< 3 > + Ai = 5L^ , -A 2 = -5L« , 



or 

T (3) = 5L (2) mod A, T (2) = -5L (1) mod A, 

where we say F = G mod A if F — G is a functional linear combination of derivatives of A. Thus 
the superintegrable system admits a nondegenerate potential V if and only if L^> = = 0. 
This last condition exactly characterizes the spaces classified by Koenigs: flat space, the 2- 
sphere, the 4 Darboux spaces, and the family we call Koenigs spaces. Thus it is a consequence 
of Koenigs' classification is that all spaces admitting 3 second order Killing tensors automatically 
admit a nondegenerate potential. 
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In establishing the above result we have not made use of all of the information obtainable 
from the symmetry integrability conditions (|12|) . In the identity A12 = the second derivative 
terms in A appear in the form Ai2(A22 — An) + • • • = 0. If A12 = then = = and we 
have one of the spaces found by Koenigs. Suppose A12 7^ 0. Then we have a new integrability 
condition 

V3Ai 2 ; 12 

Solving for ^22 — Mil from A = and substituting into = we obtain, after a straightforward 
computation, a condition of the form 

SWfJL! + S (2) /X2 + Sfl = 0, 

where 

S« = 5L (1) mod A, S {2) = -5L (2) mod A, 
S = - Lf + 3A 1 L { - 1) - 3A 2 L {2) mod A. 

Thus we have the integrability condition 

5I^Vi - 5L^H2 + {L? - 4 2) + 3A^ - 3A 2 L^) l x = 0. (15) 

When a nondegenerate potential exists then the space of solutions of (|12p is 4-dimensional 
and the values of fx, fj,\, {12, /tn can be prescribed arbitrarily at any regular point. Thus 
the integrability condition ([15]) can hold only if L« = = 0. This proves that the systems 
admitting a nondegenerate potential coincides with the potential free systems found by Koenigs. 
The same argument goes through for the case when a 2-parameter potential exists. Then the 
space of solutions of (fT2j) is 3-dimensional and the values of fx, fj,\, [12 can be prescribed arbitrarily 
at any regular point, so the integrability condition (|15|) can hold only if = = 0. This 
proves that any 2-parameter potential must be a restriction of a nondegenerate potential, a fact 
proved by a different method in pp. 



6 1-parameter potentials 

The theory of 1-parameter potentials is more complicated than that for nondegenerate and 
2-parameter potentials, due to the possible occurrence of systems with either 4 or 3 linearly 
independent second order symmetries. Suppose V(x, y) is a 1-parameter potential satisfying 
canonical equations 

V l = B 1 (x,y)V 2 , V 2 2-V 11 = B 22 (x,y)V 2 , V 12 = B 12 (x, y)V 2 . 

The integrability conditions for these equations are 

B 12(1 - {B 1 ) 2 ) - B\ - B X B\ - B l B 22 = 0, 

B\ 2 - Bf - - B\B 12 - B l B\ 2 = 0. (16) 

In special coordinates B 22 , B 11 are given in terms of B 1 by relations (jlOH . so the first equa- 
tion (fT6l) becomes 

pi 1 Dlnl / \ . \ _ \ / 1 \ 

(17) 
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unless B l = 0, in which case it becomes B 12 = Ai = 0. Note also the special case B l = ±i, 
which implies A 2 = A\B l . 

The second equation (fT6|) becomes 

A(An + \ 22 )B l + AiA 2 (l + (B 1 ) 2 ) - (A 2 + X^B 1 + XX\B\ 

+ AA 2 (B\ - 2BIB 1 ) + A 2 (B^ + B\Ai + 3B l A u - 3A 12 ) = 0. (18) 

Substituting (|17p into the expression (|10p for B 22 we find 

fi22 = (» + »g 1 )((B 1 ) 2 -l)-(g 2 1 + g 1 g 1 1 ) j (1Q) 

Now we write down the involutory system of equations that determine the second order 
symmetries of a 1-parameter potential system 

11 Ai n A 2 12 
a, =--a --a , 

2 X A v ; A 



fa 11 - a 22 



) i = l f-B 22 + 2^B' - 2^) a 12 + 2B\ 

(a n -a 22 ) 2 = -26, 
«i 2 = b, 

a? = I (-B 22 + 2^ - 2^) a 12 + fl 1 *, 

1 1 11 1 a / Ai 12 A 2 22 A 
6r = --a 2 -^^ T a + y a j, 

6 2 = -2S 1 6i - 2£}& + f i(2^ B 1 - 2^ - J3 



a iZ . (20) 



There is one additional condition, obtained by differentiating one of the Killing equations ([Tf 
that we have not made use of in obtaining the involutory system: 



= 2a 2 2 + a 22 +(4V 1 + 4V 2 ) . (21) 



A A /2 

When the indicated differentiations and substitutions are carried out, the right hand side of each 
of these equations can be expressed in terms of the variables a 11 ,a 11 -a 22 ,a 12 ,b = a\ 2 alone, 
although the expanded terms are lengthy. We think of B , B 22 as explicit given functions. 

Note that though this system is in involution, the system without the added variable b = a\ 2 
is not. This demonstrates that a symmetry is uniquely determined by the values of a , a 22 , 
a 12 and a\ 2 at a regular point; the values of a , a 22 , a 12 may not suffice. By assumption, 
the system has 3 functionally independent second order symmetries. However, the involutory 
system indicates that there may, in fact, be 4 linearly independent second order symmetries, 
but obeying a functional dependence relation. 

Now we require that the system (120p . (1211) admit 4 linearly independent solutions. Then 
at a regular point there exists a unique solution with any prescribed values of a 11 , a 22 , a 12 , b 
and the integrability conditions for the system are satisfied identically in these variables. To 
investigate the properties of this system we expand the condition (|2ip in terms of the basic 
4 variables. The result takes the form D l {x, y)a 12 + D 2 (x, y)b = 0, where the do not depend 
on the basic variables. Indeed, 

D 2 = \B X B 22 - S{B l B\ + B\) - X2B 1 + 9^ + ^{B 1 ) 2 - ^{B 1 ) 3 , 

A A A 
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with a similar but more complicated formula for D . Since the variables can be prescribed 
arbitrarily at a regular point (x,y), the only way for condition ((2TJ) to hold is for D 1 = D 2 = 0. 
We solve for B 1 B\ + B\ from the equation D 2 = and substitute this result into (|19p to obtain 
an updated expression for B 22 that is independent of A 12 : 



_ 1 -\ 2 B l + X^B 1 )' + 2,\ 2 {B 1 ) -3Ai 



B 22 {x,y) :=- 



B l \ 



(Here we are assuming B 1 ^ 0. This special case will be treated separately.) We substitute 
this updated expression for B 22 in all of the previous equations, and eliminate B\ from all 
expressions, including (|18p . The condition D 1 = is now satisfied identically with the updated 
expressions. The only remaining constraints are the integrability conditions for the symmetry 
equations (|20j) . These conditions are satisfied identically except for d\b 2 = d 2 b\ which takes the 
form E 1 (x,y)a 12 + E 2 (x,y)b = 0, where the E^ do not depend on the basic variables. Thus 
E 3 (x,y) = for the case where 4 linearly independent symmetries exist. 

In the case where the space of symmetries is strictly 3-dimensional, the integrability conditions 
will no longer be satisfied identically, since there is a linear condition satisfying by the variab- 
les a 11 , — a 22 , a 12 , b Thus, for example the condition D 1 (x,y)a 12 + D 2 (x,y)b = should now be 
considered as a constraint relating a 12 and b. 

These integrability equations for the 4-dimensional and 3-dimensional cases are rather compli- 
cated and their geometric significance is not clear, so we will pass to a simpler, Stackel transform 
approach, while making use of the partial results we have obtained via the direct integrability 
condition attack. 



7 The Stackel transform for 1 parameter potentials 
with 4 linearly independent symmetries 

To shed more light on this case we follow the approach of Section [5j That is, we choose special 
coordinates and restrict our attention to the symmetries for which a 12 ^ 0, essentially a two- 
dimensional vector space. Since B 1 and a 12 are invariant under Stackel transformations and the 
equations for the symmetries and the metric A depend only on these variables, these equations 
are identical for all metrics \jl describing systems Stackel equivalent to the original one. The 
basic equations are the symmetry conditions 

/ii 2 = 0, A = fi 2 2- Mil - 3/iiAi + 3fi 2 A 2 - (A u + A\ - A 22 -A 2 2 )^ = 0, 
A = An + A 22 + A\ + A\ = 0, 

where A = In a 12 , and the integrability conditions for the 1-parameter potential V, where 
V\ = B l V 2 . Writing B 1 = B for short and using the partial results obtained in the preceding 
section for simplification, we can write the first integrability condition (|17p for the potential as 

BA\ - A 2 B 2 + BBi , . 

Mi- Bfi 2 + D^L = 0, D = W B l + 1 + B2 + 1 > ( 22 ) 

and the second integrability condition (|18p as 

fi llf i(B 3 + 3B) - 2fi 22f iB 3 - f4(B 3 + 3B) + ^i 2 2 (3B 3 + B 5 ) 

+ mn(B 2 B l - 3#i) - 2^ 2 B 2 B U = 0. (23) 

These equations have a different interpretation than those of the last section. First of all they 
hold for 2 distinct functions a 12 whose ratio is nonconstant. Secondly, the space of solutions \x of 
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this system is 2-dimensional. Thus at any regular point (x, y) there is a unique solution /j,(x, y) 
taking on prescribed values fx, M2 at the point. Use of (|22|) and differentiation yields the linear 
expressions 



Mi = Bfi 2 - Dn, /in = (B x - DB)fi 2 + (D 2 - D x )n, 



M22 



D 2 D-B 2 

-I* + —fT^^ (24) 



B 



B 



for mi j Mill M22 in terms of \x and \x 2 . 

There are additional integrability and compatibility conditions for the system (|22p . (|22j) that 
constrain A and B. The only nontrivial integrability condition for the subsystem /ii2 = 0, 
Mi — -BM2 + Dfi = 0, is 



B 



DD 2 
B 



M + 



D 2 + 



D-B 
B 



M2 = 0. 



Since this must hold for all solutions \x we have the 2 conditions 



Pi 
B 



DD 2 
B 



Do + 



D-B 2 
B 



0. 



(25) 



The requirement that the subsystem be compatible with A = is 



El-D 2 + D l + ZA X D - 



i.e., 



£2 
B 



D 2 + D 1 + 3A X D - C = 0, 



D-B 2 
B 



D-B 2 
B 



+ DB-Bx- 3^iS + 3^2 M2 = 0, 



+ DB-B X - 2,A X B + 3^2 = 0, 



(26) 



where C = An + A 2 — A 22 — A 2 . We already knew the second of conditions (|26p but, from the 
analysis of the previous section, the requirement of a 4-dimensional space of symmetries yields 
the identity B 2 + BB\ = B(A 2 — BA X ), so D has the alternate expression 



D 



B 2 + BB X 
' B 2 + l 



This is the full set of integrability conditions. 

The key to understanding these systems is the first order condition (|22|) . We will show that 
this equation, together with other integrability conditions, implies that each system in the family 
admits a Killing vector, which is also a first order symmetry of the system. To see this, consider 
the condition that the form X = £pi + r]p 2 be a Killing vector. This condition is simply that 
the Poisson bracket of X and the free Hamiltonian TCq vanish, i.e., 



pj+p 2 

M 



0. 



Thus the coefficients of p 2 , p 2 , p\p 2 in the resulting expression must vanish: 

2£iM + £Mi + W2 = 0, 2t?2M + £Mi + ¥2 = 0, J71 + 62 = 0. 
We find that 

m = -6> £1 = m, where 2£im + f Mi + ¥2 = 0. 



In order that (|22l) be interpretable as the Killing vector requirement there must exist an integra- 
ting factor Q such that Q[i\ — QBfi 2 + QD/i = and Q = £, rj = —QB, £1 = QD/2. Further we 
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must require QD/2 = n 2 = —Q 2 B - QB 2 = £1 = Qi, and Q 2 = & = ~Vi = QiB + QB\. Thus 
we obtain the system 

(lnQ) 1 = D/2, (\nQ) 2 = B l + BD/2 

whose integrability condition is 

D 2 = 2B n +B 1 D + BD 1 , or (l + B 2 ) (B n + B 22 ) - 2B[B\ + B|) = 0. (27) 

This condition is a consequence of (|23p and the second integrability condition (|25p . Indeed, 
using d2U) to express /tin, fi\, fi 22 in terms of fi 2 and /i in (f23j) we find that the resulting 
expression takes the form F(x,y)fi 2 = 0, so F(x,y) = 0. Solving for Bi 2 in each of F(x,y) = 
and (I25p . and equating the results, we get exactly the desired condition (|27|) . 

We have shown that X = £pi + r\p 2 = Qp\ — QBp 2 is a Killing vector. Moreover X is a first 
order symmetry, since 

{X, H} = {X, Ho} + {X, V} = -Q{Vx - BV 2 ) = 0. 

Now we check the special cases B = 0, ±i. The cases B = ±i are essentially the same. 
Choosing B = —i, by interchanging x and y if necessary, we see that the second integrability 
condition for the potential gives the condition A ln(A) = 0, i.e., the condition that A is a flat space 
metric. Solving this equation, and identifying solutions that are equivalent under Euclidean 
transformations and dilations, we find three cases: 

I : A = 1, 11^ = 6^ III: X = x 2 + y 2 . 

For the first case a straightforward computation yields the Killing vector X = p 2 — ip\ and the 
potential V = a(y — ix). This is the superintegrable system [E4] in [3]. For case II the Killing 
vector is X = e~^ y+lx ^ 2 (p 2 — ip 2 ) and the potential is V = ae~^ y ~ tx ^ (corresponding to the 
superintegrable system [2214] in [1]). For case III the Killing vector is p 2 — ip\ and the potential 
is V = l/(y — ix) (corresponding to the superintegrable system [E13] in [4]). In all these cases, 
Ai 2 = 0. If B = then Ai = and the Killing vector is p±. Thus V = V{y) and all metrics \i 
Stackel equivalent to A will satisfy p\ = 0, so the fundamental equations are 

m = 0, A = A = fi22 + 3/U2A 2 - (A u +A\- A 22 - A 2 2 )\i = 0. 

These equations must have a 2-dimensional vector space of solutions, so that jjl and fi 2 can be 
prescribed arbitrarily at a regular point. Since Ai = 3fi 2 A\ 2 — Cifi = where C = An + A 2 — 
A 22 — A 2 , this implies A\ 2 = C\ = 0. Thus the spaces are contained in our earlier classification. 

Theorem 1. If a 2D superintegrable system with a 1-parameter potential admits 4 linearly 
independent second order symmetries, then it also admits a Killing vector. One of the second 
order symmetries is the square of the Killing vector. 

Now we return to the generic case where B ^ 0,±i. We show that B(x,y) = B 1 always 
factors, so that V\ = X{x)Y{y)V 2 . 

Lemma 1. 

B 12 = ^l so (JnB)i2 = and B = X(x)Y(y). 

Proof. We solve for B\ 2 = B\ 2 from (|23|) . i.e., from F(x,y) = 0. Then we solve for Bn 
from ([2?]) and substitute this result into our expression for B\ 2 . We obtain B\ 2 = B\B 2 j B . ■ 
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Since X = Qp\ — QBp2 is a Killing vector, X 2 is a second order symmetry with —a 12 = Q 2 B. 
Thus A = In a 12 = 21nQ + InB where (lnQ)i = D/2, (lnQ) 2 = B x + BD/2. From this it is 
straightforward to compute the derivatives of A in terms of B and its derivatives. We have 

B 2 + BBi B ± . „ (B 2 + BB 1 )B B 2 

A 1 = -2-^ - + — , M = 2Bi - 2 v 2 „ ^— + — , 

£? 2 + 1 B B 2 + l B 

2B22 + B 2 Bi + B\2 4(i?2 + BBi)BB 2 B\ 2 B — B 2 B\ 

A-12 = . -. 1 /^n . -xo h 



B 2 + l (B 2 + l) 2 B 2 

with analogous expressions for An 2 an d ^4i22- Next we substitute B = X(x)Y(y) into each 
of these expressions, and in the identity (|27p . and then compute = A\\ 2 — A\ 2 A\ and 
L (2) = A122 - A U A 2 in terms of X(x), Y(y). Solving for X"(x) from the identity d27]) and 
substituting this expression into L«, L( 2 ) , we find 

L« = L^ 2 ) = 0. 

Theorem 2. // 2L> superintegrable system with a 1-parameter potential aV admits 4 linearly 
independent second order symmetries then there exists a superintegrable system with nondege- 
nerate potential V(a, (3, 7) such that the restriction V(a, 0, 0) = aV and the restricted second 
order symmetries of the nondegenerate system agree with a three-dimensional subspace of the 
second order symmetries for the 1-parameter potential. 

For future use we note that the two-dimensional space of nonzero symmetries a 12 (excluding 
the zero function) does not necessarily have the property that a 12 satisfies the Liouville equation 
A(lna 12 ) = ca 12 . However, in the generic case B ^ 0,±i we have shown that there is a Killing 
vector X = ^p\ + rjp2 and that the associated a 12 = £1] from the symmetry X 2 is nonzero and 
does satisfy the Liouville equation. 

8 The Stackel transform for 1-parameter potentials 
with exactly three linearly independent symmetries 

Again we follow the approach of Section [5] and restrict our attention to the symmetries for 
which a 12 7^ 0, now a one-dimensional vector space. Since B 1 and a 12 are invariant under 
Stackel transformations and the equations for the symmetries and the metric A depend only 
on these variables, these equations are identical for all metrics [i describing systems Stackel 
equivalent to the original one. The basic equations are the symmetry conditions 

/xi 2 = 0, A = n 22 - nn - 3/iiAi + 3/i 2 A 2 - (An + A 2 - A 22 -A 2 )fi = 0, 
A = An + A 22 + ^1 + ^2 = 0, 

where A = In a 12 , and the integrability conditions for the 1-parameter potential V = V(x,y): 
V! = B 1 (x,y)V 2 , V 2 2-V 11 = B 22 (x,y)V 2 , V 12 = B 12 (x, y)V 2 . 

The integrability conditions for these equations are 
B 12 {\ - (B 1 ) 2 ) - B\ - B l B\ - B 1 B 22 = 0, 

B\ 2 - Bf - B\ x - B\B 12 - B 1 B\ 2 = 0. (28) 

In a special coordinate system £> 22 , B 11 are given in terms of B 1 by relations (jlOp . so the first 
equation (i28|) becomes 

, Bl +B l B\ - W l (A 2 - B l A l ) , , 

Mi-BV2 + £>/z = 0, D = -± \ + {B i)2 -■ ( 29 ) 
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(The cases B 1 = 0, ±i cannot occur because we have seen that they lead to 4 independent 
symmetries.) The second equation becomes 

H(HU + V22)B l + W/ u 2 (l + (B 1 ) 2 ) - (pi + LL 2 2 )B l + nmB\ 

+ ihh{B\ - 2B\B l ) + ^(B^ + B\A X + 3B 1 A U - 3A 12 ) = 0. 

Substituting (I29p into the expression (|10h for B 22 we find 



B 



22 



} + fB^f-l)-{Bl + B-B\) 



B 1 



and, as before, 



B 



12 



A key equation for this approach is the integrability condition (|15p of Section [5l derived from 
consideration of A12 = 0: 



SL^Vi - 5L^^ + - 4 2) + 3AiL« - = 0. 



(30) 



Recall that LW = A X \ 2 - A X2 A X , = A i22 -A i2 A 2 . Here £W = = is the condition that 
the system admits a nondegenerate potential. We have already determined all such 1-parameter 
systems: they are just the restrictions of nondegenerate potential systems to a single parameter. 
Our interest is in finding systems that are not simply restrictions, if such systems exist. Thus 
we require that the coefficients of /zi, (jL 2 , fj, in ([30]) are all nonvanishing. (Since (f30|) must admit 
a 2-dimensional solution space, if any one coefficient vanishes then all vanish.) Setting B = B l 
for short, we rewrite conditions (I29p and (|30f) in the form 



Ml 



/'i 



Bfi 2 + 

L (2) 



B 2 + BB 1 -3B(A 2 -BA 1 ] 



B 2 + l 

(2) 



■n = o, 



L(i) 



A*2 + 



L W_42 J+ 3^ iL(1 )_3^ 2L (2) 
5L(i) 



"A 4 



0. 



(31) 



Since there are solutions with arbitrarily chosen values of fi 2 , \x at a point, these equations must 
be identical: /ii — S//2 + -D/z = 0, where 



B 



£ 2 + RBi - 3B(A 2 - BA X ) _ LY> - L 2 ' + 3A 1 L^ - 3A 2 L^ 

WTi ~ 5lw 

Use of /ii2 = and differentiation yields the linear expressions 



H! = BfjL 2 - Dfi, /in = (Bi - DB)n 2 + (D 2 - D^fi, 



H22 



D. 



D 2 

~B^ + 



(32) 
(33) 



D-B 2 



B 



for fj,i, fin, fj, 22 in terms of (J, and il 2 . Just as in the preceding section, there are additional 
integrability and compatibility conditions for this system that constrain A and B. The only 
nontrivial integrability condition for the subsystem \i\ 2 = 0, (ii — Bli 2 + Dll = 0, is 



B 



DD 2 
B 



M + 



D 2 + 



D-B 2 
B 



1*2 = 0. 
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Since this must hold for all solutions fi we have the 2 conditions 

D 2 \ DD 2 [D-B 2 \ 

0, D 2 + = 0. 



B J 1 B ' \ B 

The requirement that the subsystem be compatible with A = is 
D 2 „\ f D — B 



i.e., 



D z + D 1 + 3A 1 D-C)fi+ I — j^ + DB-B 1 -3A 1 B + 3A 2 ) m = 0, 



§ - D 2 + Di + 3AiL> - C = 0, - ^ +DB-B X - 3A Y B + 3A 2 = 0, 
B B 

where C = An + A\ — A 22 — A\. Note that the second condition is already satisfied. 

Since we no longer assume a 4-dimensional space of symmetries, we can no longer deduce 
directly that D = —2(B 2 +BiB)/(l+B 2 ). However we shall see that the assumption L^'L^ 2 ' / 
and identities (|32|) . (f33j) imply this result. 

Indeed we can substitute B = L^/L^ into all identities and express them in terms of 
A = In A 12 and its derivatives. The verification that 

2 , - 4 2) + 3A 1 LW - 3A 2 L^ 
-2(B 2 + B 1 B)/(1 + B 2 ) = -± 2 ^ 2 = D (34) 

is cumbersome, but straightforward. One uses the identities ([33]) and A = (and its derivatives) 
to express Ai 222 , Am 2 , A i22 , in terms of derivatives of A of lower order (or in some cases equal 
order with the number of y derivatives < the number of x derivatives) and substitutes into both 
sides of the desired identity until equality is evident. Indeed the only fourth order derivatives 
occurring in (|34|) are A\ 222 and A\\\ 2 . We can solve for these derivatives in terms of strictly 
lower order derivatives of A from the identities (|33|) and A12 = 0. Then we substitute these 
expressions back into the left and right sides of (|34p to obtain lower order expressions. Similarly 
we can solve for A\ 22 from Ai = and A 22 from A = and substitute back. After several 
iterations the identity is verified. 

Now we are in exactly the same situation as for the 4-dimensional symmetry case in the 
previous section. Since D = —2(B 2 + B\B)/{1 + B 2 ) the system must admit a Killing vector 
X = Qpi — QBp 2 . The integrability condition for Q is D 2 = 2Bn + B\D + BD\ which 
follows from pip and other identities listed above. Thus Lemma Q] is valid for this case, so 
B = L^/L^ = X(x)Y(y), i.e., B factors, and the proof of Theorem [2] is valid. Thus = 
= 0, which is a contradiction. Thus the only 1-parameter superintegrable systems with 
exactly three second order linearly independent symmetries are those for which = = 0. 

Theorem 3. If 2D superintegrable system with a 1-parameter potential aV admits exactly 3 li- 
nearly independent second order symmetries then there exists a superintegrable system with non- 
degenerate potential V(a, (3, 7) such that the restriction V(a, 0, 0) = aV and the restricted second 
order symmetries of the nondegenerate system agree with a three-dimensional subspace of the 
second order symmetries for the 1-parameter potential. 

9 Structure of systems with 4 linearly independent 
second order symmetries 



We begin by proving a fundamental duality result for 2D systems with 4 linearly independent 
second order symmetries and a 1-parameter potential. (In fact, however, it is not hard to show 
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that every system with 4 independent second order symmetries must admit a Killing vector and 
a 1-parameter potential.) We have worked out the fundamental structure equations for these 
systems in Section [6j In the generic case they are 

A*i2 = 0, a 12 {pn - /i 22 ) + 3ma\ 2 - 3[i 2 al 2 + (a\ 2 - a^)n = 0, a\f + o} 2 \ = °> ( 35 ) 

Bo + BBi B(al 2 - Ba} 2 ) 

Hi - Bfi 2 + Dfi = 0, D = -2-^ -i = -2- 



5 2 + l a 12 (5 2 + l) ' 

where V\ = BV2 an d B satisfies the integrability conditions 

(B 2 + 1)(B 11 + B 22 )- 2B(B 2 + B 2 ) = 0, BB i2 = BiB 2 . (36) 

Using one of the alternate expressions for D, we can write the first equation (|35p in the more 
obviously dual form 

(B 2 + 1) (//i - Bfi 2 )a 12 + 25 (Sal 2 - 4 2 )^ = 0. (37) 

(In the form ([35]) . (|36|). ([37]) these equations apply not only to characterize the generic case 
but also the special cases B = 0,—i.) The interpretation of these equations is that /i = A ^ 
is a metric of a superintegrable system corresponding to the potential function B provided A 
and B satisfy equations (j35j) . (|36|) . (I37j) for some 2-dimensional space of harmonic functions a 12 . 
Then the solution space of all simultaneous solutions fi is also 2-dimensional and the solutions 
give precisely the metrics of systems Stackel equivalent to A. 
The duality result is the following: 

Theorem 4. If fi(x,y), a 12 (x,y), B(x,y) satisfy fl35]), (J36J), dHTJ) then 

~, s 12 f x + iy ~ ix -y\ -12/ \ (x + iy -ix-y 

H{x, y) = a — 7=-, , a (x,y) = 1 



B(x,y) = i 



V2 ' V2 J' K , * J r V 72 ' v/2 

D / x+iy —ix—y 



also satisfy these conditions. Thus, the roles of metric and symmetry can be interchanged if 
the potential function undergoes an appropriate Mobius transformation. A second interchange 
returns the system to its original state. 

Proof. This is a simple consequence of the chain rule and the relations 

a l 2 = 72^ 1_ ^ 2 ^' fl 2 2 = 73(^1 - ^2), 

with analogous relations between B and -B. Under the duality the first condition ()36[) obeyed 
by i? maps to the second condition (|36p obeyed by -B, and vice- versa. The theorem also applies 
to the special cases where B = or B = —i. Under the duality the special cases switch roles. ■ 

Theorem 5. Every 1-parameter superintegrable 2D system with 4 linearly independent second 
order symmetries is Stackel equivalent to a superintegrable system on a constant curvature space. 

Proof. We consider the generic case B ^ 0,—i first. Then every such 1-parameter superin- 
tegrable 2D system with 4 linearly independent second order symmetries and metric \jl = A 
corresponds to the system of equations (f35|) . (i36j) . (|37|) where the a 12 range over a 2-dimensional 
space and there is a 2-dimensional space of solutions [i, corresponding to the metrics of systems 
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Stackel equivalent to the [i = A system. Moreover there is a Killing vector X = £pi + r\p 2 
such that (Zq 2 = £77 from the second order symmetry X 2 is nonzero and satisfies the Liouville 
equation. It follows that the dual metric £ = aj 2 is of constant curvature. The dual system also 
describes a Stackel equivalence class of 1-parameter sup erintegr able 2D systems with 4 linearly 
independent second order symmetries, so there must exist a corresponding Killing vector y and 
nonzero function 6q 2 that is harmonic and satisfies the Liouville equation. Since the duality 
maps all solutions fi of the original system one to one onto symmetries b 12 of the dual system, 
there must exist a unique solution fj, = v of the original system such that v = b^ 2 . Since 6j 2 
satisfies the Liouville equation v is the metric of a constant curvature space. This means that 
the system with metric A is Stackel equivalent to the constant curvature system with metric v. 

Now we consider the special cases. If B = — i then all metrics are flat, so the statement 
of the theorem again holds. If B = then X = p\ so the associated Oq 2 = 0. In this case, for 
any nonzero harmonic a 12 from the 2-dimensional space of symmetries of the original system we 
have that a 12 is a metric for the dual system with potential function B = —i. Now such metrics 
are flat, and as we have shown in the proof of Theorem Q] there is then a Killing vector X such 
that the symmetry 6q 2 from X 2 is nonzero and satisfies the Liouville equation with c = 0. Since 
the duality map is one to one and onto there must exist a metric v in the original system such 
that v = 6q 2 . This means that v is a flat space metric. ■ 

We continue to require that the space of second order symmetries is of dimension 4 and now 
investigate the space of third order constants of the motion for this system: 

2 

K. = a kjl (x 1 ,x 2 )pkPjPi + b e (x 1 ,x 2 )pi, 

k,j,i=l 

which must satisfy {7Y,/C} = 0. Here a k i l is symmetric in the indices k, j, i. 

Theorem 6. For a 2D superintegrable system with 1-parameter potential and a 4- dimensional 
space of second order constants of the motion, the dimension of the space of strictly 3rd order 
constants of the motion is at most 4. 

Proof. We give the essence of the proof, leaving out some of the details. The conditions on JC 
are 

2af = -3((lnA)ia™ + (lnA^a^), i + j, 

3af 1 + of = 3 - ((In A)*o iij + (In A)^'), i / j, 

2(a\ 22 + af 2 ) = -(InAHa 122 - (lnA) ia m - (lnA) 2 a 222 - (lnA) 2 a 112 , 

2 22 
b\ + b 2 = ^\a s21 V s , bj = lj2a s ^V s -^J2(ln\) s b s , j = l,2, (38) 



2 ^ 2 

s=l s=l s=l 



and 



]T 6-l/ s = 0. (39) 

s=l 

The a fcjl is just a third order Killing tensor. The b e must depend on the parameter in the 
potential V linearly, so we have 

2 

b (xi, x 2 ) = ^2 f j (x 1 ,x 2 )V j (xi,x 2 ). 



2D Superintegrable Systems with 1-Parameter Potentials 



19 



Using the standard notation V\ = B 1 V 2 , V 12 = B 12 V 2 , V 11 = B n V 2 , V 22 = B 22 V 2 with B 1 = B, 
we find from that Bf 1 = f, f 2 = -Bf 
Further 



bi = fiV 2 + fV 12 , b L 2 = f 2 V 2 + fV 22 , 

b\ = -(Bf)^ - BfV 12 , b\ = -{Bf) 2 V 2 - BfV 22 . 

Thus 



b\ = (f 1 + fB 12 )V 2: b\ = {f 2 + f 22 )V 2 , 

b\ = -{Bft + (B 1 + BB 12 )f)V 2 , b\ = -{Bf 2 + (B 2 + BB 22 )f)V 2 , 

Now we assume B / and (1 + 2A)(1 - 2 A + 4A 2 ) / 0. (The last inequality can always 
be achieved by a rescaling of the coordinates x\, x 2 , if necessary. We will treat the special 
case B = separately.) Substituting these results into the defining equations (|38p and solving 
for fi, f 2 and a 122 , we obtain 

fx = ^a m + \a™ + .-., (40) 

3BX iin 3A 222 , 

f 2 = -a -a H , (41) 

Jl 1 + 2A 1 + 2A ' y 1 

a 122 = J^_ a 112 _ 1 ^ fl 222 (42) 

1 + 2A B{l + 2\) ' v ; 

where the omitted terms are linear in /. Substituting these results into the remaining 5 condi- 
tions, we can recast them in the form 



h 

.in 



(43) 

a l 2 Ll = -3s + • • • , (44) 

af 2 = ■ ■ ■ , (46) 
oi 12 = «, (47) 

where the omitted terms are linear in /, a 111 , a 222 , a 112 . (Here we have introduced a new 
variable s = a} 12 -) The system (|40p -(|48 p is not in involution because we have not expressed 
the derivatives s\, s 2 in terms of /, a 111 , a 222 , a 122 , s. Differentiating dHJ with respect to xi, 
using the fact that (a2 U )i = ( a } 11 )2, arid differentiating ([13"]) with respect to x 2 , we can solve 
for si- Similarly, differentiating (|45|) with respect to x 2 , using the fact that (of 22 ) 2 = (a| 22 )i, 
and differentiating (|46|) with respect to x±, we can solve for s 2 . Thus we obtain 



«! = •••, (49) 
s 2 = -3s + • • • , (50) 

where the omitted terms are linear in /, a 111 , a 222 , a 112 . We conclude that the full system 
with 5 functions s, /, a 111 , a 222 , a 112 is in involution. Since the space of Killing vectors is 
one-dimensional we can always add a multiple aX to any 3rd order symmetry (where a is the 
parameter in the potential) and obtain another 3rd order symmetry. Thus, by adding such 
a multiple, we can prescribe / = 0, s, a 111 , a 222 , a 112 at a regular point and there is at most one 
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strictly 3rd order symmetry taking on these values at this point. Hence the space of 3rd order 
symmetries is at most 4-dimensional. 

In the special case B = we find b 1 = /V2, b 2 = Ai = 0, so a 222 = 0. A similar argument 
to the preceding shows that the system with 5 functions s, f, a 111 , a 112 , a 122 is in involution. 
Again the space of strictly 3rd order symmetries is at most 4-dimensional. ■ 

Theorem [6] gives important information about the structure of the Poisson algebra generated 
by the first and second order symmetries of this superintegrable system As we know, a basis 
for the second order symmetries can be taken in the form Cq = 7i, C\, £2, X 2 , where X is 
the Killing vector. Clearly the set of 4 third order symmetries CqX, C\X, C2X, X s is linearly 
independent. Hence it must be a basis. 

Corollary 1. For a 2D superintegrable system with 1-parameter potential and a 4-dimensional 
space of second order constants of the motion, the dimension of the space of 3rd order constants 
of the motion is exactly 4. 

From this result we see that the algebra generated by Co, C\, C2, X is closed under the 
Poisson bracket operation. Indeed the Poisson brackets {X,jCj} are second order symmetries, 
so we have 

{X,£,j} = ojoCq + aj\C\ + 0^2^.2 + CLj^X 2 , 

where the ajk are constants. The Poisson brackets {Ci,Cj} are third order symmetries, so by 
Corollary [T] they can be expanded in the form 

{£i,£j} = bijfiCoX + bij^\C\X + bij^^X + bij^X 3 . 

Additional commutators and relations between the constants follow easily from the Jacobi iden- 
tity. 

This is not the full story for these Poisson algebras, because the four linearly independent 
generators Cq, C\, C2, X 2 must be functionally dependent. We need to understand the form of 
this dependence. The key to this is the space of 4th order symmetries. 

Next we investigate the space of fourth order constants of the motion. Here a constant of the 
motion 

2 2 
T= a ekjt {x,y,z)p e p k p j pi+ ^ b mg (x, y, z)p mPq + W{x, y, z), 

l,k,j,i=l m,q=l 

must satisfy {7i,J-} = 0. Again a £k ^ 1 , b mq are symmetric in all indices. 
The conditions are 
2 

a p = _ 2 ^y™(lnA) s , (51) 

2 

4af 11 + af = -eJ2 « S% '( ln *) s , * ^ h (52) 

s=l 

2 2 

3ap + 2af j = - £ a siii (ln A) s - 3 ]T a s ™(ln A) s , i + j, (53) 



S=l 8=1 

2 2 



2b? + bf = 6A a^Vs ~ hS3 ( ln A )*> » ^ 3, (54) 

s=l s=l 
3 2 

bf = 2A ^ a siii V s - bSj (^ A) s , 



s=l s=l 
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and 

2 

Wi = X^2b si V s . (55) 

8=1 

Note that the a ^ is a fourth order Killing tensor. We require the potential V to be superinte- 
grable 1-parameter, and that the highest order terms, the in the constant of the motion, 
be independent of the parameter in V. The b mq must depend linearly and W quadratically on 
this parameter. We ignore the arbitrary additive constant in W. 
We set 

Then, listing only the first derivative terms in the explicitly, conditions (|54p become 

2/i 12 + / 2 n = -- - , (56) 

2/ 2 12 + A 22 = • • • , (57) 

fi 1 = ■■■ , (58) 

f¥ = ---- (59) 

From the integrability condition d X2 W\ = d xl W2 for equations ([55]) we obtain the condition 

BXifl 1 + A2A 12 - BX^l 2 - \ x f 22 = • • • (60) 

which does not involve the variables a l i ki . From conditions (|51|) - (|53p we obtain 

a\ U1 = ---, (61) 

a? 22 = ■■■, (62) 

Aa\ 112 + al ul = --. , (63) 

44 222 + af 22 = ... , (64) 

3a\ 122 + 2a 1 2 ll2 = ... , (65) 

34 122 + 2a} 222 = ... , (66) 

where we list only the first derivative terms in the a l ^ ke explicitly. 

Theorem 7. For a 2D superintegrable system with 1-parameter potential and a A- dimensional 
space of second order constants of the motion, the dimension of the space of strictly 4th order 
constants of the motion is at most 9. 

Proof. We give the basic ideas, leaving out some of the details. The proof is by parameter 
counting. There are 8 variables, 

f ll f 12 f 22 1111 1112 1122 1222 2222 

j , j , j , a , u , a , a , a 

The number of possible first derivatives of these variables is 16, and they are subject to the 
11 conditions ()56|) — (|66[) . The number of possible second derivatives of the variables is 24, and 
they are subject to the 22 conditions, obtained by differentiating the original conditions with 
respect to x and to y, Similarly, the number of possible third derivatives of the variables is 32, 
and they are subject to the 33 conditions, obtained by taking all second derivatives of the 
original conditions. At this point the system is in involution, since all third derivatives of the 
variables can be expressed in terms of lower order derivatives. The total number of variables is 
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8+16 + 24+32 = 80 and the number of conditions is 11 + 22 + 33 = 66, so there are 14 parameters 
remaining. However, there is one more third order condition than third derivatives, so we can use 
the 33 conditions to eliminate all of the third derivative terms and obtain a condition involving 
at most second derivatives. Explicitly, this relation is 

d u & + A 2 (-Bd 12 + ^d 22 \ (EHD + Ai (d 12 - jd n } {MD 

- ydiaPD + ^d 12 m = 0, 

where we have indicated the identities to be differentiated by their equation numbers. The 
resulting condition is 

6-BAA2 (-Ball 11 + ^ a 22 U ) + 56AAi f-f a?f 2 + a 22 2 22 ^j + 3£A(-A 2 + BX^a 1122 

+ AA 2 ( - Ba 1112 + all 12 ) + AAi( — a 1222 + Ba 1222 ) = • • • , (67) 

where we have listed explicitly only the second derivative terms in the a-variables, not the f- 
variables. It is straightforward to verify that the second order part of condition (|67p is nonzero 
and that it is not obtainable by taking linear combinations of first derivatives of conditions 
(|6ip - (|66p . Thus it is an independent second order condition. Differentiating this condition with 
respect to x we get a third order condition with leading terms 

6-BAA2 (-Ba 1111 + \a\^\ + 56AA! (-fog? 2 + a 222 A + 3BX(-X 2 + BX 1 )a 1122 



+ AA 2 ( - Ba 1112 + a\l 12 ) + AAi ( - ajff + Bal 



12 ) 



which is not obtainable by taking linear combinations of second derivatives of conditions (|6ip - 
(|66p . Thus the maximum number of parameters for this system is reduced to 13. However we 
can add a 2 times any linear combination of a basis of four 2nd order symmetries to a 4th order 
symmetry (where a is the parameter for the potential V). Thus the maximum dimension for 
the space of strictly 4th order symmetries is 9. ■ 

As we pointed out earlier, the basis of 4 second order symmetries Cq, C\, C 2 , £3 = X 2 
must be functionally dependent. This dependence must occur at the fourth order. Indeed we 
can form 10 products '■ j < k}, from these symmetries, all strictly fourth order. Since 

the maximum dimension of the space of fourth order symmetries is < 9 it follows that the set 
{CjCk} is linearly dependent. There can only be one such linear dependency relation, and we 
call it the Casimir. This relation, and the second order and third order structure relations given 
previously, completely determine the structure of the Poisson algebra. 



10 Discussion and outlook 



We have filled in a gap in the literature by clarifying the structures of the quadratic algebras 
corresponding to 2D second order superintegrable systems with 1-parameter potentials. For such 
a system the number of linearly independent second order symmetries is either 3 or 4. If it is 3 
the system is simply a restriction from a nondegenerate potential system. It is 4 if and only if 
it admits a first order symmetry. For 4 second order symmetries the quadratic algebra is closed 
at orders three and four from dimensional considerations: The number of terms obtainable from 
the basis first and second order symmetries exceeds the maximal dimension of the space of third 
order and of fourth order symmetries. Every such system is Stackel equivalent to a system on 
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a constant curvature space. In 3D and higher dimensions the relation between systems corre- 
sponding to degenerate and nondegenerate potentials is much more complicated [25\ |2"o] |2"7] but 
the tools developed here should be applicable. (Indeed in 3D the potentials can depend on up to 
4 parameters, neglecting the additive constant. Nondegenerate (4-parameter) potentials always 
have 6 linearly independent second order constants of the motion although the number of func- 
tionally independent constants is 5. The functional relation between the constants is of order 8 in 
the momenta, Not all 3-parameter potentials are restrictions of nondegenerate potentials. Some 
such systems admit a closed symmetry algebra with some generators of order 4 [26]. There 
are up to 8 second constants of the motion for 2-parameter potentials.) It is also important to 
extend the structure analysis to superintegrable systems of order greater than two. 

Finally the O-parameter (free space) case deserves more attention, even though it was treated 
by Koenigs [19]. If we consider a true O-parameter potential system as one in which it is not 
possible to admit a nonconstant potential, then it follows from our results and those of Koenigs 
that no such system exists. Koenigs' proof used complex variable techniques, very different 
from the methods used here. Further the last stages of his argument where he shows that no 
superintegrable systems exist other than those he has found are somewhat murky. Koenigs does 
not consider potentials or the Stackel transform, and as we have shown, these are very important 
in the theory. Thus there are good reasons for a new approach to uncovering the structure of 
these systems. 
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